Terahertz surface plasmons in optically pumped graphene structures 
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We analyze the surface plasmons (SPs) propagating along the optically pumped single-graphene 
layer (SGL) and multiple-graphene layer (MGL) structures. It is shown that at sufficiently strong 
optical pumping when the real part of dynamic conductivity of SGL and MGL structures becomes 
negative in the terahertz (THz) range of frequencies due to the interband population inversion, the 
damping of the THz SPs can give way to their amplification. This effect can be used in graphene- 
based THz lasers and other devices. Due to relatively small SP group velocity, the absolute value 
of their absorption coefficient (SP gain) can be large, substantially exceeding that of the optically 
pumped structures with the dielectric waveguide. The comparison of the SGL and MGL structures 
shows that to maximize the SP gain the number of GL layers should be properly choosen. 
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I. INTRODUCTION 

Optical excitation of graphene can result in the inter- 
band population inversion [3, 01 ■ At sufficiently strong 
excitation, the interband emission of photons can pre- 
vail over the intraband (Drude) absorption. In this case, 
the real part of the dynamic conductivity Re be- 
comes negative. Due to the gapless energy spectrum of 
graphene [3| , Re can be negative at relatively low fre- 
quencies w, for instance, those in the terahertz (THz) 
range. This effect can be used in graphene-based THz 
optically pumped lasers with the Fabri-Perot resonators 
and the resonators based on dielectric or slot-line waveg- 
uides 4-6]. As was previously pointed out [l| and ana- 
lyzed p} , apart from the lasing associated with the stim- 
ulated generation of electromagnetic modes, the stimu- 
lated generation of different plasmons (with their conver- 
sion into electromagnetic radiation) can also be of prac- 
tical interest. The plasma excitations in graphene were 
considered, in particular, in Refs. [8l-[l3j. In this paper, 
we study the propagation and amplification of surface 
plasmons (SPs) in graphene structures under nonequilib- 
rium conditions when Re a u < due to optical pump- 
ing. We consider a single-graphene layer (SGL) and 
a multiple-graphene layer (MGL) on a substrate. As 
demonstrated recently (see the review paper |14| and the 
references therein), the MGLs, which constitute stacks 
of disoriented (non-Bernal stacked) GLs, exhibit similar 
electron and optical properties as individual GLs, while 
the electron momentum relaxation time in such MGLs 
can be extremely long. The latter circumstance implies 
that the intraband absorption of photons and plasmons 
can be effectively diminished. This makes the MGLs 
very prospective for different optoelectronic devices con- 
sidering that the net dynamic conductivity along GLs 
increases with increasing their number in the MGL struc- 
ture. The MGL structures can be particularly effective 
in lasers in which the THz modes are supported by the 
external resonator or waveguide THz photodetec- 

tors [l^HL7t . and transit-time THz oscillators [H, [l9j . 
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FIG. 1: Schematic view of the SGL/MGL structure under 
consideration. 



However, SGL and MGL structures can play the dual 
role: the imaginary parts of their dynamic conductivity 
provides the mode localization near the SGL or MGL 
(i.e., the formation of SPs), while the real part provide 
absorption or amplification of plasmons. Due to this, the 
MGL structures as SP waveguides are not always supe- 
rior over the SGL structures. Since the thickness, d, of 
the MGL structure with even relatively large number, K, 
of GLs is rather small in comparison with the wavelength 
of the plasmons under consideration, we shall consider 
the SGL and MGL structure alike. The real MGL usually 
include a highly conducting bottom GL (at the interface 
between the substrate and other GLs) with rather high 
electron density [3, H(J ■ Although this bottom GL can 
be in some way removed, we will take it into account. 

It is assumed that the GL or MGL structure is illumi- 
nated from the top by light with the energy of photons 
Ml. Under optical excitation, the electron and hole densi- 
ties exceed substantially their equilibrium values. Due to 
this, one can consider the electron and hole systems under 
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FIG. 2: Frequency dependence of the real part of propagation index Re (p) (solid lines) and the SP group velocity normalized 
by the speed of light (dotted lines) in SGL structures at different temperatures (T = 300 K and T — 77 K) and different 
quasi-Fermi energies (1 -s% = meV, 2 -ej = 10 meV, and 3 -£p = 20 meV) for SGL structures with different substrate 
refraction index (n = 1.0 and n = 3.4). Dashed and dashed-doted lines in figures for n = 3.4 correspond to the ranges where 
the effects of spatial dispersion might be essential. 



consideration as characterized by the quasi-Fermi ener- 
gies ±£p , respectively, and the effective temperature T. 
A distinction in the Fermi energies in GLs with different 
indices is due to the attenuation of incident optical radi- 
ation associated with its absorption in the GLs closer to 
the structure top: Ep < — s F . The electron- hole 
system in the bottom GL in MGLs is characterized by 
the unified Fermi energy ep determined by the doping 
of this GL. If the characteristic time, To, of the emission 
of the optical phonon by an electron or a hole is much 
shorter than the time of pair collisions, the photoexcited 
electrons and holes emit a cascade of optical phonons 
and occupy low energy states in the conduction and va- 
lence bands, respectively. In this case, the contribution 
of optical excitation to the heating of the electron-hole 
system is small, so that the effective temperature T is 
close to the lattice temperature T\ [2l|. In the opposite 
case, the photoexcited electrons and holes are thermal- 
ized, and their effective temperature is determined by the 
pumping photon energy Ml and the rate of electron and 
hole energy relaxation. In such a situation, the effective 
temperature can be elevated, so that somewhat stronger 



optical pumping might be needed to fullfil the condition 
Re a„ < [12. 



II. DYNAMIC CONDUCTIVITY OF SGL AND 
MGL STRUCTURES 

The net dynamic conductivity in the lateral direction 
of a structure with K GLs at the signal frequency u> can 
be presented as the sum of the contributions of the in- 
dividual GLs dif' (k — 1,2, ...,K) and the bottom GL 



°» = Y.oM+o*. (1) 

k=l 

Considering the expressions for crffl and obtained pre- 
viously (see, for instance, Refs, [HHHl), one can arrive 
at the following: 
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Here, e is the electron charge, rg and r are the electron 
and hole momentum relaxation times in the bottom and 
other GLs, respectively, ft is the reduced Planck constant, 
kg is the Boltzmann constant, and 



G( £ ,e') 
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(4) 



In the case of SGL structure, one should put K = 1 and 
cr^ = 0. For the MGL structures without the bottom 
GL, one can use Eq. (3) formally setting a% = 0. 

(k) 

The quasi-Fermi energies in the GLs with k > 1, e F , 
are mainly determined by the electron (hole) density in 
this layer Ej., with e F oc V® (at sufficiently strong 
degeneracy of the electron and hole systems), i.e., by the 
rate of photogeneration Gq ' by the optical radiation at 
the k— th GL plane. Considering the attenuation of the 
optical pumping radiation due to its absorption in each 
GL, one can obtain that e F can be expressed via the 

quasi-Fermi energy in the topmost GL ep" = Epf , which, 
in turn, is a function of the the intensity of incident 
pumping radiation Iq. Considering the MGL structures, 
the e F versus sp* dependence was found as previously Q . 



III. DISPERSION EQUATION FOR SP. 

Using the the Maxwell equations, in which the current 
associated with electron-hole system is given by of 

ix = 0, j v =0, jz = a u> E z S(x), (5) 

The components of the electric and magnetic fields, E = 
(E x ,0,E z ) and H = (0,H y ,0) can be presented in the 
following form: 



E x oc E z oc H y oc exp 
at x > 0, and 



[ i—y/l — p 2 x+i— pz— iuit ) , (6) 
V c c J 



E x oc E z oc H y oc exp I —i — \J n 2 — p 2 x + i—pz — iuit 



at x < 0. Here, uj is the signal (plasmon) frequency, c 
is the speed of light, n is the refraction index (generally 
complex) of the substrate, p is the complex propagation 
index (which relates to the wavenumber q z = pui/c), the 
axis x is directed perpendicular to the structure plane, 
the axis z corresponds to the in-plane direction along the 
SP propagation, and S(x) is the Dirac delta function. The 
latter describes the localization of the plasmon near the 
SGL or MGL structure plane. Considering Eqs. (6) and 
(7), from the Maxwell equations we obtain the following 
dispersion equation for SPs: 
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(8) 

Equation (8) can be derived from that obtained previ- 
ously using a different approach [TJ|, HH. What is im- 
portant for the present consideration is that a u is deter- 
mined by both nonequilibrium electron and hole compo- 
nents with, in particular, Re (cr^) < 0. Solving Eq. (8) 
with respect to p, one can find the propagation index 
Re (p) and the SP absorption (amplification) coefficient 
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FIG. 3: Dispersion of SPs in a SGL structure at T = 300 K 
(solid lines) and T = 77 K (dashed lines) with the same as in 
Fig. 2 values of the quasi-Fermi energy. 
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FIG. 4: Frequency dependences of SP absorption coefficient 
21m (q) in a SGL structure at T — 300 K (upper panel) and 
T = 77 K (lowe panel) and different quasi- Fermi energies (1 - 
e£ = meV, 2 -e| = 10 meV, and 3 -e T F = 20 meV). Markers 
show the dependences for equilibrium electron-hole systems 
calculated using the derived formulas as well as using [121 ] for 
intrinsic and doped SGL structures. 



21m (q) — 2Im(puj/c) or the SP gain -21m (q). In the 
cases when n = 1 (SGL or MGL structure levitating in 
the free space), Eq. (8) obviously yields 



IV. RESULTS OF CALCULATIONS 
A. SP spectra and SP damping/amplification 

Figure 2 presents the frequency dependences of the real 
part of the propagation index Re (p) oc Re(q) shown by 
solid lines calculated for SGL structures with different 
substrate refraction indices n at different temperatures T 
and different values of the quasi- Fermi energies ep (i.e., 
for different pumping intensities). Dotted lines in Fig. 2 
correspond to the frequency dependences of the group 
velocity of SPs normalized by the speed of light in vac- 



FIG. 5: Spatial distributions of the real parts of electric 
and magnetic fields in SGL structures at different frequencies 
(upper panel) and different substrate refraction indices (lower 
panel): T = 300 K. e T F = 20 meV, and r = 10 ps. 



uum c. In our model we disregarded the effect of spa- 
tial dispersion on the dynamic conductivity of SGL and 
MGL structures, i.e., the dependence of <r u on the SP 
wavenumber q z in Eqs. (1) - (4). This is valid if p < c/vw 
or uj < Re(q z )vw [III; where vw — 10 8 cm/s is the char- 
acteristic velocity of electrons in holes in graphene. In the 
case of SGL structures with n = 1, these inequalities are 
satisfied in the frequency range under consideration (see 
Figs. 2 and 3). However, as follows from Fig. 2 at larger 
n, the validity of the dependences obtained is limited by 
relatively low frequencies (u/2ir < 5 — 8 THz for n = 3.4 
at T = 300 K). One can see from Fig. 2 that the SP 
group velocity can be negative (at T = 77 K) at elevated 
frequencies. This corresponds to backward waves [24l[25j 
and occurs at the frequencies at which the imaginary part 
of the dynamic conductivity is determened primarily by 
the interband transitions. Figure 3 shows the spectra 
of SPs, i.e., uj versus Ke(q z ) dependences in the SGL 
structures at different temperatures at the same pump- 
ing conditions as in Fig. 2. 

Figure 4 shows the SP absorption coefficient 21m (q z ) 
as a function of frequency calculated for different tern- 
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FIG. 6: Frequency dependences of SP absorption for SGL 
structures with different substrate refraction indices n: 1 - 
n = 1.0, 2-n = 3.4, 3-n = 3.4+i0.01, and 4-ra = 3.4+i0.05, 
and 5 - n = 3.4 + iO.l (T = 300 K, e F = 20 meV). 



peratures T and different values of the quasi-Fermi en- 
ergies Bp. One can see that an increase in Ep leads to 
widening of the frequency range where Im (q z ) < and a 
marked increase in the absolute value of the absorption 
coefficient in this range. The dependences indicated by 
markers were calculated for equilibrium SGLs using the 
above formulas which in the absence of optical pum pin g 
at n = 1 coincide with those obtained previously [12j . 
The markers correspond to intrinsic (with the Fermi en- 
ergy ep = 0) and doped SGL structures. The depen- 
dences for doped GL structures at the equilibrium cor- 
respond to the Fermi energies ep in doped SGLs at the 
equilibrium equal to the quasi-Fermi energies £p in the 
optically pumped SGLs. Clear distinctions in these de- 
pendences in the equilibrium and pumped GLs are asso- 
ciated with the contributions of both electrons and holes 
to the negative dynamic conductivity due to the inter- 
band transitions (in the latter). 

Figure 5 shows the spatial distributions (in the direc- 
tion perpendicular to the SGL plane) of the SP electric 
and magnetic fields in the SGL structures with different 
substrate refraction indices. 



B. Role of substrate and intraband absorption 

Figure 6 shows the frequency dependences of the SP 
absorption coefficient calculated for SGLs with the sub- 
strates with different real and imaginary parts of the re- 
fraction index (different absorption in the substrate) . As 
follows from Fig. 6, the contribution of the substrate to 
the SP absorption can be insignificant at the realistic 
values of the imaginary part of the substrate refraction 
index. In particular, in the case of the substrate made of 
undoped Si (Im (n) ~ 3 x 10~ 4 [111), the imaginary part 
of the refraction index can be smaller than those used in 
the calculations of curves 3- 5 in Fig. 6. 

Figure 7 shows the frequency dependences of the real 
part of the propagation index and the SP absorption co- 
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FIG. 7: Real part of SP propagation index (upper panel) 
and absorption coefficient (lower panel) in SGL structures 
as functions of frequency: 1 — e F = meV, r = 10 ps , 
2 - e T F = meV, r = 0.1 ps, 3 - e T F = 20 meV, r = 10 ps, 
and 4 - ep = 20 meV, r = 0.1 ps. 



efficient. These dependences were calculated for SGL 
structures with different electron and hole momentum 
relaxation times (t = 10 ps and r = 0.1 ps) for differ- 
ent £p, i.e., for different pumping intensities assuming 
that T = 300 K. As seen from Fig. 7, in the SGL struc- 
ture with a relatively long momentum relaxation time 
(t = 10 ps), the SP absorption coefficient changes its sign 
at moderate values of e'p ( ep < 20 meV) . However when 
t = 0.1 ps, the SP absorption coefficient does not change 
the sign at least at Ep ~ 20 meV (see below), although 
it is markedly smaller than in the equilibrium conditions 
(without pumping). However, at elevated pumping in- 
tensities (elevated values of the quasi-Fermi energy) the 
absorption coefficient in the SLG structures even with 
short momentum relaxation times can become negative 
although at higher frequencies (see Fig. 8). 

C. Comparison of SGLs and MGLs 

In Figs. 9 and 10, we compare the real parts of the 
SP propagation index, absorption coefficients, the SP 
group velocities in the optically pumped SGLs and MGLs 
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FIG. 8: SP absorption coefficient as a function of the quasi- 
Fermi energy in SGL structures with short momentum relax- 
ation time (r = 0.1 ps) at 1 — lo/2it — 15 THz, 2 — uj/2-k = 
20 THz, 3 - uj/2-k = 25 THz. 



FIG. 10: Frequency dependences of SP absorption coefficient 
at £§• = 20 meV for SGL structure with 1- r = 10 ps, 2- 
r = 1.0 ps, and 3- r = 0.54 ps and for MGL structure with 
K — 20 (dashed line) and t — 10 ps. 
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FIG. 9: Frequency dependences of real part of SP prop- 
agation index (solid lines) and normalized group velocity 
(dashed lines) for SGL structure (K = 1) and MGL struc- 
ture (K = 20) at e T F = 20 meV. 



FIG. 11: Spatial distribution of real parts of the SP electric 
and magnetic fields in a SGL structure and MGL structures 
with K = 20 at T = 300 K, Sp = 20 meV, r = 10 ps, and 
cj/2tt = 6 THz. 



(e£ = 20 meV). As seen, the SP characteristics in SGLs 
and MGLs are markedly different. In particular, SGLs 
exhibit substantially stronger (about 50 times in the 
peaks) amplification (higher plasmon gain) in a wide fre- 
quencies range (from 4 to 8 THz) at t = 10 ps. This is 
attributed to markedly different values of the SP group 
velocities in SGLs (relatively small plasmon group ve- 
locity) and in MGLs (in which this velocity is relatively 
high), as shown in Fig. 9, which, in turn, is due to differ- 
ent net electron and hole densities and different widths 
of the SP electric field localization that is clearly seen in 
Fig. 11. An example of the dependence of the SP absorp- 
tion coefficient in the structures with optically pumped 
GLs on the number of the latter is demonstrated in 
Fig. 12. Thus, at the same values of the optical pump- 
ing intensity and the electron and hole momentum relax- 
ation time, SPs in SGL structures can exhibit stronger 



amplification than MGL structures. However, as demon- 
strated [IH, the momentum relaxation time in epitaxi- 
ally grown MGLs can be rather long, so MGL structures 
might be preferable in applications. Taking into account 
the real possiblity of a large difference in the electron 
and hole relaxation times in SGL and MGL structures, 
in Fig. 10, we demonstrate that the SP amplification in 
the MGL structure with K = 20 and t = 10 ps can ex- 
ceed that in the SGL structure with smaller r. Hence, to 
achieve the SP maximum gain, the number of GL layers 
should be optimal. 

A decrease in the SP gain with increasing K is in con- 
trast to the behavior the electromagnetic modes in opti- 
cally pumped MGL structures with the dielectric waveg- 
uide considered by us previously [f| . For comparison, the 
pertinent dependence taken from Ref. @ for the waveg- 
uide thickness L = 5 fim. (at the same quality of GLs and 
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FIG. 12: Absorption coefficients of SPs (solid line) and an 
electromagnetic mode (dashed line) in a dielectric waveguide 
(DW) as functions of the number of GLs K. 



pumping conditions) is shown in Fig. 12 as well. 



tion in GLs and the substrate. It was demonstrated that 
at sufficiently strong but realistic optical pumping, the 
SP absorption coefficient can be negative (so that the 
gain is positive) in the certain range of frequencies. The 
absolute value of the SP absorption coefficient (plasmon 
gain) in SGL structure can be fairly large and markedly 
exceed the gain of electromagnetic modes in the dielectric 
waveguides with optically pumped SGL or MGL struc- 
tures. In contrast to the cases of amplification of electro- 
magnetic modes in the structures with optically pumped 
SGL and MGL and dielectric waveguide, the SP amplifi- 
cation weakens with increasing number of GLs K. How- 
ever, due to possibility of rather long momentum relax- 
ation times of electrons and holes in MGL structures, 
the SP maximum gain can be achieved at optimal values 
of the nimber of GLs. The effect of SP amplification in 
optically pumped SGL and MGL structures can be used 
in THz lasers (with the conversion of SPs to the output 
electromagnetic modes) . The variation of the SP charac- 
teristics by optical pumping might be useful in different 
GL-based plasma- wave devices. 



V. CONCLUSIONS 

We studied the SPs in optically pumped SGL and MGL 
structures. Using the developed model, we calculated 
the SP dispersion relations, spatial distributions of their 
electric and magnetic fields, and frequency dependences 
of the absorption coefficient as functions of the optical 
pumping and factors determining the intraband absorp- 
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